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Abstract 



Relativistic corrections to ttie hyperfine splitting are calculated for the triplet 
2^ Si state of the helium isotope He^. High precision variational wave functions 
are employed, where the electron-electron correlations are well accounted for. 
Due to the unknown nuclear structure a comparison with the experimental 
result is performed via He^"*" hyperfine splitting. Surprisingly, the second 
order contribution due to Fermi interaction, in spite of the additional small 
factor m/M, is very significant. 
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The calculation of higher order relativistic and QED effects in few electron systems 
is a long standing problem. Various measurements of transition frequencies have reached 
the precision of few ppb, while no theoretical predictions are yet so accurate The 
usual approaches, which incorporate most of relativistic effects from the beginning, MCDF 
for example, are not capable to include electron-electron correlations in a complete way, 
thus their predictions are not sufficiently accurate. The alternative approach starts from 
the Schrodinger equation and incorporates relativistic effects perturbatively in the effective 
hamiltonian. The principal advantage is the simplicity and high accuracy of nonrelativistic 
wave functions, which allows for precise calculations of higher order relativistic effects. The 
cost one pays in the perturbative approach is the complexity and singularity of the effective 
hamiltonian. Historically, the helium fine structure splitting was the first application of 
the higher order effective (Breit-like) hamiltonian. Douglass and Kroll in with the help 
of Bethe-Salpeter formalism derived effective operators at order ma^, which together with 
the second order correction from leading Breit terms gave a complete ma^ contribution 
to helium fine structure of 2Pj levels. It was the first disprove of the common believe, 
that Breit Hamiltonian could not be used at higher orders of perturbation calculus. 20 
years later Khriplovich and coworkers essentially rederived their result, without being 
aware of this former work, using the time ordered perturbation technique. With these 
effective operators they calculated fine structure splitting in positronium. The extension 
of the perturbative approach to S-states is a highly nontrivial task, see for example 0. It 
is due to high singularity of effective operators, which matrix elements between S-states 
diverge, so one may question about the correctness of the perturbative approach. In spite of 
these divergences, it is possible to implement the perturbative approach to S-states as well. 
The idea lies in the regularization of Coulomb interaction, in such a way that all matrix 
elements become finite. All contribution coming from high-momenta are incorporated as 
local interaction, proportional to Dirac delta or it's derivatives. The coefficients are found 
by matching for example, the scattering amplitude and they depend on the regularization 
parameter. However, this dependence cancels out when all terms contributing to energy at 
specified order are summed up. As a test, one may rederive |Q Dirac energy levels at order 
ma^. This approach was successfully applied first in the calculation of positronium hfs [Q. 
This result was confirmed independently by 3 groups |^ , two of them were applying directly 
Bethe-Salpeter equation, and this really serves as an independent check of this method. 
Little later, one performed calculations of higher order relativistic corrections to helium 
2^Si level. The fact, that it is a triplet state, was a great simplification, since the wave 
function vanishes at ri = r2. Very recently, Yelkhovsky |^ has derived complete set of 
operators for the ground I^Sq state. The numerical calculation of their matrix elements 
is still under way. In this work we study He^ hyperfine splitting, which was measured 
very precisely by Rosner and Pipkin in p|. He^ hfs is a nice example of difficulties in the 
bound state QED and perturbative approaches. The Fermi interaction, as given by Dirac 
S^{r) function is already quite singular. Incorporation of further relativistic effects leads to 
even more singular operators, which have to be properly handled. In the work ||^ we have 
shown the practical way to handle singular effective operators for helium atom. Coulomb 
interaction is smoothed at small distances with the use of some parameter A, to avoid small r 
nonintegrable singularities. The key point of this approach lies in the fact, that dependence 
on A cancel out between all matrix elements, and correct energy levels are restored in the 
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limit A oo. High electron momenta does not contribute here, because it is again triplet 
S-state, so one does not need to construct counter terms for matching scattering amplitude. 
Surprizingly, we have found a very large second order correction due to Fermi interaction. 
This correction requires a separate treatment and will be dealt with at the end of this work. 

Hyperfine splitting in He^ of 2^5*1 state is due to the interaction of electron and helion 
(nucleus) magnetic moments. In the 2^5*1 state both electron spins are parallel and sum to 
S* = 1 in contrast to the ground state were S = 0. Magnetic moment of helion comes mainly 
from the neutron particle. It means that helion g-factor is negative and hyperfine sublevels 
are inverted with respect to hydrogen, the upper one has S + I = 1/2 and the lower one 
S + I = 3/2. Therefore, by a hyperfine splitting one means here i^hfs = E{l/2) — E{3/2). 
According to the perturbative approach the general expression for the hyperfine splitting up 
to the order m is: 
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where a, k are anomalous magnetic moments of the electron and the helion respectively. 
However, it is not commonly accepted such a notion for a nucleus. The relation of k with 
the magnetic moment of the nucleus with charge Z e is n = 2 (l + A;) Z e/(2 M) I. Masses m 
and M are of the electron and helion respectively. The expectation values of H^^^ and H^^^ 
vanish in 2^5*1 state, but they contribute in the second order, last term in Eq. (p. H^^^ is 
a Breit hamiltonian in the nonrecoil limit: 
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where r = ri — r2 and r = \r\. Any possible recoil corrections at order ma^ are included 
in E^ec- In spite of the small additional factor m/M their contribution is significant. This 
effect could be associated to the mixing of 2^51 and 2^5*1 states caused by the Fermi 
interaction, (-ffhfs) is delta-like term with the coefficient given by the two-photon forward 
scattering amplitude. It is the same like in hydrogen and strongly depends on the nuclear 
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structure. It also automatically includes nuclear recoil effects and inelastic contribution 
(nuclear polarizability). (i^Ms) has been considered in detail for the case of hydrogen and 
muonic hydrogen. However, there are no sufficient experimental data available for helion, 
He^ nucleus, therefore we were unable to estimate these contributions. Moreover, it would 
be incorrect to apply this correction for point-like nucleus, because high energy photon 
momenta are involved, where nucleus, could not be approximated as a point like particle. 
Therefore we will leave this contribution unevaluated, and at the end for the comparison 
with an experiment, will subtract the appropriately scaled hydrogenic value for hyperfine 
splitting The last term H^^l includes spin dependent operators, which contribute at order 
ma^. It is only this term, which derivation was not performed so far in the literature and 
is presented here. The detailed description on derivation of effective hamiltonian, reader 
may find in former works, for example in [Q. There are four time ordered diagrams, which 
contributes to if^fs ^^"^ presented in Fig 1. The first two are the same as in hydrogen, 
other two are essentially three body terms. One derives the following expressions, which 
corresponds to these diagrams. 
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It is worth noting that Vi ans V2 are the same as in hydrogen. Matrix element of if^fs 
last term in Eq. (|l]) are separately divergent at small ri and r2. However these divergences 
cancel out in the sum. It is because the high electron momentum contribution is absent in 
the nonrecoil limit. For hydrogen this sum is equal to the expectation value of 7 A on Dirac 
wave function. We introduce now the following regulator A to the electron-nucleus Coulomb 
interaction 
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in all hamiltonians in Eq.(|l|), as well as in the nonrelativistic one. This leads to the following 
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Once the interaction is regularized, one can calculate all matrix elements and take the 
limit A —>■ 00. As a first step, using formulas from [Q, we rederived the known relativistic 
correction to hfs in hydrogen 
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where n is a principal quantum number. It agrees with that, obtained directly from the 
Dirac equation. Since for helium, all matrix elements could be calculated only numerically, 
we will transform effective operators to the regular form, where A could be taken to infinity 
before the numerical calculations. The initial expression for a complete set of relativistic 
corrections in atomic units is (with implicit A regularization and excluding recoil): 
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where Sn = (-f^hfi)) we used the following formulas for hfs of ^5*1 states: 

{an ■ CTi) = (cr„ ■ (cTi + (T2)/2) = -3 , 

{al ai < {a, + a^)'' ) = -2 g^-'' , 

for symmetric and traceless Q*-'- There is also a one loop radiative correction 
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which is similar to that in hydrogen. It will not contribute to the special difference between 
the helium and hydrogen-like helium hfs, therefore we will not consider it any further. The 
initial expression is rewritten to the regular form, where A regularization is not necessary. 
The operators in second order terms £a are transformed with the use of 
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This transformation leads to new form for £'j^ and £'^, such that 
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In the numerical calculations of these matrix elements we follow the approach developed by 
Korobov 0. The S wave function is expanded in the sum of pure exponentials 

N 

= E^* (e-"'^'^-^'^^^-^''^ - in ^ r2)) , (33) 
1=1 

with randomly chosen ai,(3i,''yi in some specified limits. This basis set has been proven to 
give excellent results for the nonrelativistic energy and the wave function. Moreover, its 
simplicity allows for the calculations of relativistic corrections. With basis set = 1200 we 
obtained the nonrelativistic energy (without the mass polarization term piP2 f^/M) 

E = -2.1752293782367913057(1) , (34) 

slightly below the previous result in [l^. Expectation values of Dirac delta function without 
and with the mass polarization term are correspondingly: 

( 4 TT {5^{ri) + 5^{r2))) = 33.184142630(1) , (35) 
( 4 TT (^^(ri) + 5^(r2)))MP = 33.184152589(1) . (36) 

The last one gives the leading hfs in helium, which is 

mM 



Ems = 2 Za^ |1 + fc| (1 + a) (tt (S'in) + 5%r2)))MP ~ 6 740 451kHz, (37) 



where we use values of physical constants from Ref. |TT[ with one exception Numerical 
results for Sx with X = A,B,D,N are presented in Table I. The inversion of if — in 
Sa is performed in the similar basis set as for 2^Si wave function, however the nonlinear 
parameters had have to be properly chosen, to obtain a sufficiently accurate result. Namely, 
if < X, y, Z < 1 are independent pseudo random numbers with homogeneous distribution, 
then: 

a = A2 + , (38) 
(3 = {B2-Bi)Y + B,, (39) 
7 = (C2 - Ci) Z + Ci. (40) 
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Parameters A, B,C and n are found, by minimization of the second order term with regu- 
larized Dirac delta on both sides. The inversion of H — E in Sb is performed in the basis 
set of the form 
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Unfortunately, we have not been able to get a reliable number for Sd with this numerical 
approach. The reason is that operators in Sd are so singular, that this basis set gives a 
very pure convergence. The result presented in Table I, is obtained analytically within 1/Z 
approximation, namely we neglected completely electron-electron interaction and corrected 
this value by factor (7r(5^(ri) -f- (5^(r2)))/9. The estimated uncertainty is of the order of 10%. 
The total contribution of ma^ term (without recoil -Erec) to hfs is 

^(6) = |i + fc| -^a^ 200.9844(5) = 2171.440(5) kHz , (42) 

what could compared to the leading Fermi contact interaction in Eq. (|37D , E^^^'/Ey^fs ~ 
0.000322. Between all the contributions to He^(2^S'i) hyperfine splitting in Eq.(|ID, H^^\ 
essentially the nuclear structure contribution, requires input from the nuclear physics to be 
reliable evaluated. This is the reason we do not present final theoretical predictions for hfs, 
to compare with the precise measurement in [§] 

EMs(He) = 6 739701.177(16) kHz . (43) 

Instead, we can compare our result indirectly by subtracting the ground state hfs of helium 



ion as measured in 13 



Ehfs(He+) = 8 665 649.867(10) kHz , (44) 
by composing the following difference 

AEe., = i?Ms(He) - ^ ^(^^(^i) + ^'(^^)))mp ^^^^(He+) = -38.998(19) kHz . (45) 
4 o 

In this way, nuclear structure contribution, of order ma^ cancels out, as well as the leading 
Fermi contact interaction. What remain are electron-electron correlation effects. Theoreti- 
cal predictions for this difference are (see Table I) 

AErei = \l + k\ 1.8795(5) = 20.306(5) kHz. (46) 

We do not associate here the uncertainty due to higher order terms. A disagreement of 
the theoretical result in Eq. (|46D with the experimental one in Eq. ( ^51) indicates that the 
presented calculation is incomplete. However the old calculations of relativistic corrections 
in [14] using an approximate DHF wave function, led to the result which is in an agreement 
with experiment AEom = —32(22) kHz. This agreement is only illusive, since there is a 
correction discovered by Sternheim |15|, which is denoted here by -Erec- It is the second 
order contribution due to Fermi interaction H^^^ in Eq. (|^) 
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Erec = ( i^Ms _ ff)/ ^hfs ) • (47) 

It could be understood as a recoil correction since it includes additional small factor m/M. 
Sternheim noticed that H^^^ mixes 2^5*1 and 2^5*1 states, what together with small energy 
difference between these states leads to strong enhancement of this recoil correction. Stern- 
heim result which includes only 2^5"! intermediate state is 5E = —66.7(3) kHz, what nicely 
would explain the discrepancy of 59.3 kHz. However, the inclusion of higher excited states, 
will lead to the infinite result. Moreover, this correction is partially included in {hI^I), and 
only after proper subtraction it becomes finite. Hopefully, all divergence or in other words 
dependence on regularization parameter A cancel out in this particular difference between 
helium atom and helium ion hyperfine splitting. Therefore instead of Eq. ( ^T]) we calculate 
the difference 
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Further evaluations proceed in the similar way as for Sa_- With the help of Eq.(^) one sub- 
tracts algebraically all divergences and the remaining finite terms are calculated numerically. 
The result is 

AE,,, = -- ( ^-^^^-^y aHAA91.0(l) = -\1 + k\ — 5.5971 = -60.471 kHz (49) 
2\ 3mM J ^ ^ ' ' M ^ ' 

The complete theoretical predictions are the sum of Eqs. (EBI) and 



AEt^ = -\l + k\ 3.71756 = -40.165kHz (50) 

m M 



It nicely agrees with the experimental value -38.998 kHz from Eq. (^). The small difference 
of 1.167 kHz could be associated to higher order QED corrections. In summary we have 
calculated relativistic and a dominating recoil correction to helium hyperfine splitting. Due 
to unknown nuclear structure we considered such a difference of hfs of helium atom and 
He+ ion, that cancels out this nuclear structure contribution. The recoil correction, namely 
the second order term in the Fermi interaction give a significant contribution to hfs, due to 
mixing of singlet and triplet S-states. Was this correction included, a good agreement with 
experiment is achieved. 
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FIGURES 




FIG. 1. Time ordered diagrams contributing to helium hyperfine structure at order ma^. 
Dashed line is a Coulomb photon, the wavy line is the transverse photon, the thicker vertical 
line denotes nucleus, two other electrons. 
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TABLES 
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TABLE L Numerical results for contributions at order m a to helium hyperfine structure. The 
factor 24 in the above comes from Breit correction Eq.([l8D with n = 1 times 3/4 from spin algebra 
times 
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